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A class of well-motivated models of inflation end by producing cosmic strings. The 
current status of efforts to calculate and observe the signals from such models are outlined, 
with a particular emphasis on cosmic strings, and on the Cosmic Microwave Background 
signal. 

§1. Inflationary cosmology and cosmic strings 

Current observations are consistent with a remarkably simple model of the very 
early universe: inflation. 1 ) At a time of about 10 -36 sec all the energy density of the 
universe is held in the form of a homogeneous scalar field, the inflaton. The potential 
energy density relaxes so slowly that it acts like an effective cosmological constant, 
causing the expansion of the universe to accelerate. In this accelerating or inflating 
phase, the quantum fluctuations in the scalar field are amplified, blown up in scale, 
and frozen in as a scale-invariant spectrum of density fluctuations. 2 ) The inflaton 
eventually approaches the minimum of its potential energy density, oscillates around 
it, and decays into other particles, signalling the start of the conventional hot big 
bang. Thanks to the period of inflation the universe is large, flat, and smooth, apart 
from the scale-free perturbations laid in by the field. 

In some models there is a phase transition in a companion set of fields towards the 
end of inflation. 3 )' 4 ) A thermal phase transition in a theory with a topologically non- 
trivial set of ground states leads to the formation of defects, 5 ) extended structures 
made from spatially varying configurations of the scalar (and possibly also gauge) 
fields. At the end of inflation the inflaton replaces the temperature as the control 
parameter, and defects form equally well. 3 )' 6 )' 7 ) Hybrid models can be naturally 
accommodated in the framework on supersymmetry and Grand Unification. 8 ) 

Of particular interest are effectively one-dimensional structures, cosmic strings, 9 )' 10 ) 
where the energy density is concentrated into a thin tube. Strings have the crucial 
property: they decay just fast enough to maintain a constant density relative to the 
rest of the contents of the universe, a property known as scaling. This property 
is not shared by all defects: in particular monopoles (point-like) and domain walls 
(two-dimensional) end up dominating the energy density by today (13 GYr) unless 
their mass scale is much lower than that of inflation (10 15 GeV). 

One-dimensional structures are also singled out as the fundamental constituents 
of matter in string theory, which allows for infinitely long and relatively light (com- 
pared with the Planck scale) "elementary" strings to be created in the hot big 
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bang. 12 ) These strings also scale, and hence the study of cosmic strings is now 
an important feature of the cosmology of elementary strings. 

Although the interactions of strings are model-dependent, all kinds of strings 
have gravitational fields, and if they are massive enough they add to the gravita- 
tional perturbations from inflation. The search for strings is therefore partly a search 
for deviations from the standard inflationary model. This work focuses on the sig- 
nals from strings, as the predictions of inflationary cosmology are now part of the 
standard lore, 133 ) with a particular emphasis on the Cosmic Microwave Background 
(CMB) signal where there is an interesting interplay between the inflation and cosmic 
string parameters. It also concentrates on "ordinary" cosmic strings: those without 
currents, without junctions, and with unit reconnection probability, to remain to a 
certain extent complementary to other recent reviews. 10 ) 

§2. Cosmic strings 

Cosmic strings are linear distributions of mass-energy in the universe, which may 
have any length. They are characterised by two numbers: the mass per unit length 
fj,, and their tension T. In their fundamental form, the tension and the mass per 
unit length are equal, fi = T, which means that there is a boost invariance along the 
length of the string and they are relativistic. 

Like particles, they may be elementary or composite - made from more funda- 
mental constituents, in this case scalar and possibly gauge fields. If they are made 
from fields, they are classical solutions of the field equations localised into tubes of 
energy density, related to solitons. 13 ) Classically, elementary strings have zero width, 
and solitonic strings have a non-zero width set by a mass scale in the classical field 
equations which we denote M cs . Quantum fluctuations smear out an elementary 
string to a size which is the inverse string scale M s _1 , but do not significantly correct 
the size of a solitonic string in a weakly coupled field theory. 14 ) 

In both cases, the classical dynamics are simple for strings whose curvature 
radius is much larger than the width (either M" 1 or M s _1 ): the acceleration is pro- 
portional to the local curvature, with certain constraints. Strings couple to other 
fields; certainly to gravity, with a strength Gfi, where G is Newton's constant, and 
fx ~ M 2 S or Mg. For solitonic strings whose mass per unit length is set by the 
Grand Unification scale (about 10 15 GeV), G/i is in the range 10~ 7 - 10" 6 , while for 
elementary strings the effective four-dimensional tension depends on the details of 
the compactification of the extra dimensions. In the benchmark KKLMMT warped 
compactification scenario, 15 ) which includes a mechanism of inflation, phenomeno- 
logical considerations put the string tension in the range 10 -10 < Gfi < 10 -6 , 16 ) 
with a similar range for unwarped compactifications. 12 ) 

Strings also couple to massive states: solitonic strings are constructed from the 
fields of particles with mass M cs , and elementary string excitations come in an infinite 
tower of states with masses M 2 ~ nM 2 , where n is an integer. One can also regard 
string loops as massive states, produced by self-intersections of the string. 

An unsolved problem in the cosmic string scenario is the eventual destination of 
the energy in the string network. It has been addressed in different ways: originally 
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by numerical simulations of classical relativistic (Nambu-Goto) strings 17 )' 20 )' 21 ) and 
then by simulations of a classical field theory with string solutions. 22 )~ 24 ) The simi- 
larities and differences between the two approaches is easily seen in Figure 1, which 
shows a horizon-sized volume from an example of each. Thanks to the scaling prop- 
erty, the snapshots could be taken at any time. The strings are divided between 
"infinite" strings, stretching across the box, and loops. The infinite strings rather 
similar in curvature radius and length density, but loops are almost entirely absent 
in the field theory simulation. The reason for this is that all the strings' energy 
density goes into massive modes of the fields, which do not show up in the figure. In 
the Nambu-Goto approximation, this channel does not exist, and instead the energy 
of the infinite strings goes into small loops. These loops are assumed to decay slowly 
into gravitational radiation. 

It should be noted in neither set of simulations has gravitational radiation been 
included, which is vital in the traditional string scenario to determining the average 
size of the loops when they are produced: gravitational radiation reaction on the 
infinite strings is supposed to set this scale. In the scenario based on the classical 
field theory, it is argued that the gravitational fields are 0(Gfi), and therefore the 
gravitational radiation reaction is negligible compared with the reaction from the 
massive modes, which has an 0(1) effect on the string energy density over a Hubble 
time. 24 ) 




Fig. 1. String networks from simulations performed using in classical field theory (left) and in 
the Nambu-Goto approximation (right). Each box represents an approximately horizon-sized 
volume at any epoch in the history of the network. Note that the "infinite" strings, stretching 
across the box, are rather similar in curvature radius and density, but that loops are almost 
entirely absent in the field theory simulation. 
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§3. Observational signals from strings 

There are many ways to search for strings in the universe, but as outlined above 
there is a source of uncertainty in the predictions of the scenario. Perhaps the most 
well-explored avenue exploits the gravitational fields of the strings, which perturb 
the matter and radiation around them. The amplitude of these perturbations is 
controlled by the dimensionless parameter 8irG[i, where G is Newton's constant and 
\i is the string mass per unit length. The factor of 8tt comes from Einstein's equations. 
Hence gravitational perturbations of amplitude 10 -6 - 10 -5 are to be expected from 
strings in Grand Unified Theory, right in a very interesting observational range. 

The best-developed calculations are those for the power spectrum of the Cosmic 
Microwave Background fluctuations. Strings generated the fluctuations both by stir- 
ring up the photon-baryon fluid in the surface of last scattering, and by perturbing 
the CMB photons on their journey towards us. 26 )' 27 ) Standard codes for computing 
the CMB temperature and polarisation power spectra 28 ) are easily adapted to ac- 
commodate sources of energy-momentum, and there is good understanding of how 
the sources depend on length scale and time, derived from the combination of 3D 
numerical simulations with modelling at various levels. The results from numerical 
simulations 30 ) -35 ) depend somewhat on the modelling, but the differences can be 
simply parametrised in a phenomenological model of string evolution. 36 )~ 41 ) They 
are often quoted in terms of /io, the fraction of the power due to strings at a mul- 
tipole t = 10, and depending on the modelling and the dataset the upper limit 95% 
confidence limit upper bound is currently in the range 2 - 10%. 41 )~ 43 ) 

The gravitational fields of the strings likewise perturb the matter in the universe, 
although strings are much more efficient at making CMB perturbations than matter 
density fluctuations. This is because strings are a source of vector and tensor grav- 
itational perturbations as well as scalar ones: all three produce CMB fluctuations, 
but only scalar perturbations affect the density. The same codes which compute the 
CMB power spectrum can also supply the power spectrum of linear density pertur- 
bations from strings, showing that a string component consistent with the CMB is 
negligible except at very small scales. 37 ) String-induced matter perturbations are 
rather different from those sourced by inflation, as they are non-linear from the out- 
set: matter is drawn in behind a moving string into a thin sheet or wake. 44 ) This 
means that standard linear perturbation theory is inaccurate. The dominance of 
the string perturbations at very small scales means that the galaxy formation may 
proceed very differently even with a small component of strings (as measured by the 
CMB power). In particular the bias (the square root of ratio of the galaxy power 
spectrum to the matter power spectrum) is subject to uncertainty, and therefore 
directly applying the galaxy linear power spectrum 45 ) to constrain the matter power 
spectrum in models with strings is not reliable. Constraints derived from observa- 
tions of the Baryon Acoustic Oscillation 46 ) in large-scale structure surveys are not 
so sensitive to the bias, and therefore should be applicable to models with strings. 

The gravitational fields of strings can also produce gravitational lensing of galax- 
ies, quasars, and other distant compact sources. 47 )~ 51 ) A gravitational lens produced 
by a straight string consists of a pair of unamplified and undistorted images, quite 
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different from those produced by galaxies and clusters of galaxies, and it would seem 
that lens searches are an excellent way to look for strings. However, string lenses 
must be separated from the background signal of pairs of similar galaxies. 48 )' 52 ) 
Furthermore, the maximum separation of the images is approximately 5(G/i/10 -6 ) 
arc sec, which is to be compared with the resolution of ground-based surveys and 
the angular size of galaxies at redshifts of z ~ 1, both of which are around 1 arc 
sec. Future radio instruments such as the Square Kilometer Array 53 ) could produce 
interesting constraints from surveys of Compact Radio Objects, 51 ) perhaps even as 
low as Gfi < 10 -9 . 

The gravitational signal with the largest theoretical uncertainty is gravitational 
radiation. 54 )~ 60 ) Gravitational radiation is emitted both from infinite strings, from 
the collision of small amplitude waves and from oscillating string loops. The power 
from a string loop of any size is ' P = rGn 2 , where r is a 0(100) shape-dependent 
factor: the total power per unit volume from a string network is therefore propor- 
tional to the number density of loops. This number density depends very strongly 
on assumptions made about the principal energy loss mechanism for the network. 
Classical field theory simulations show that there is 0(1) loop per horizon volume 
surviving for one Hubble time, 24 ) while the traditional cosmic string scenario argues 
that the number per horizon volume is much larger, (rG[i)~ 3p in a universe whose 
scale factor depends on time as a a f. The root of this difference is that in classi- 
cal field theory simulations, strings can decay into the massive particles of the field 
from which they are constructed, whereas the traditional scenario assumes that this 
channel is eventually closed by the strings becoming smooth, and that gravitational 
radiation takes over. 

Similar remarks apply to other forms of massless radiation: both solitonic and 
fundamental strings may couple to massless pseudoscalar fields, the prime example 
being axion strings. 61 ) Massless pseudoscalar radiation is also emitted by oscillating 
strings, 62 - 1 although there is disagreement over the power and the spectrum. 63 )' 64 ) 
Classical field theory simulations again show 0(1) loop per horizon volume. 65 ) 

Solitonic strings in a realistic theory incorporating the Standard Model must 
couple at some level to known particles, and so should be a source of cosmic rays, 
perhaps even Ultra-High Energy cosmic rays (UHECRs) with energies above 10 10 
GeV 66 ) if their mass scale is high enough. Again, the power in cosmic rays is uncer- 
tain, but there are essentially two scenarios: classical field theory simulations tell us 
that all the energy of the string network goes into massive particles, and therefore 
there is a potentially large fraction of the total available energy going into extremely 
energetic Standard Model particles, 22 )' 67 ) which can show up both as UHECRs and, 
via a cascade of decays and interactions, as GeV-scale 7-rays. From cosmic ray 
observations there are strong but model-dependent constraints on the energy injec- 
tion rate and hence on the mass density in cosmic strings. 68 )~ 72 ) In the traditional 
scenario, particles are produced only near cusps (rare regions of the string moving 
close to the speed of light) and consequently the constraints are much weaker. The 
decay products may still be of interest for producing baryon asymmetry through 
out-of-equilibrium decays 73 ) or as an extra source of dark matter beyond standard 
thermal freeze-out. 74 ) 
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§4. Cosmic strings from high-energy physics 

4.1. Solitonic string solutions in the Abelian Higgs model 

The simplest gauge field theory showing a string solution is the Abelian Higgs 
model. 75 ) It consists of a complex scalar <f> and a gauge field A^, with a Minkowski 
space action 

S = jd*x (d^*D^ + V(</>) + \f» v F^ , (4-1) 

where the covariant derivative is = — ieA^ and the potential energy density 
is V = |A (|0| 2 — 0q) 2 - The fi e ld equations are 

D^ + X^l 2 - 4)^ = (4-2) 
d^F^ - ie(<j)*D u <P - D U <P*<P) = 0. (4-3) 

The ground state can be written as<p = e za 4>o, A^ = 0, or gauge transformations 
thereof. There are excitations around the ground state, which correspond to a scalar 
of mass m s = V2X<pQ and a vector of mass m v = e<f>Q 

There exist cylindrically symmetric finite energy static solutions of the following 
form 

<j> = M(r)e i9 , Ai = —Mr)k, A = 0, (44) 

er 

where 6 is an azimuthal coordinate and r a radial one. By substituting the ansatz 
into the field equations, one can discover numerical solutions of the form shown in 
Fig. 2. 
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Fig. 2. Numerical solutions for the Nielson-Olesen vortex, for parameters A = e = 1. Shown is the 
dimensionless energy density p/(e(f>o), the magnitude of the scalar field /, and the azimuthal 
component of the gauge field Ag, plotted as a function of the distance from the origin in units 
of the vector mass e<f>o- 

One can see that the the energy density is confined to a region whose size is of 
order the inverse vector mass, (e^o) -1 , and that it contains a tube of magnetic flux 
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oriented in the z direction, B z = A' e {r)/er. The string mass per unit length is found 
to be 

IX = 2irB(\/e 2 )cp 2 , (4-5) 
where B is a slowly varying function of its argument, approximately given by 76 ) 

B(B)~l L ° 4/3 °' 195 ' 10 " 2 < ^f 1 ' (4-6) 
^ {P) - \ 2.4/ln(2//3), < 1(T 2 . 1 j 

4.2. Inflation models with strings 

The inflation scenario supposes that the energy density of the universe was 
dominated by a homogeneous gauge singlet scalar field s (see e.g. Ref. 133)) with 
potential V(s). The key observational quantities are the power spectra of the metric 
perturbations, broken down into scalar and tensor components V s ,t(k), and their 
tilts n s t = dYo-Vst/dlnk. In standard slow-roll inflation driven by a single field, 



*>.<*> 1 1/1 
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> v ^ k ) - A— 

aII=k ^ 2 ™\ 



(4-7) 

aH=k 



(l-6e + 2 V )\ aH=k , nt = -2e, (4-8) 



where e = m\V n '/(2V), r] = mpV" /V are the slow-roll parameters, which must be 
small during inflation. All quantities are evaluated at "horizon crossing" aH = k, 
when the wavelength of the mode in question has grown to the Hubble length. 

In hybrid inflation models the inflaton can be coupled to the string scalar 
field 4 )' 7 ) 

V(s, 0) = V 8 (s) + |A(|0| 2 - 4>l) 2 + hs 2 \4>\ 2 . (4-9) 

For large s, the local minimum of (j> is at (j> = 0, and the potential energy density is 
V s (s) + l\<$. Inflation can take place provided V s is sufficiently flat, and terminates 
at a critical value s 2 = A<^q//i, where the <f> field becomes unstable, rolls towards its 
minimum (where V = 0) and the string-forming phase transition occurs. 

When looking for flat potentials, it is natural to explore supersymmetric theo- 
ries. A supersymmetric hybrid inflation model with a U(l) gauge symmetry can be 
constructed from three chiral superfields S, <P and <P, with U(l) charges 0, 1 and — 1 
respectively, and superpotential 8 ) 

W = kS($$-M 2 ). (4-10) 

Adopting lower case letters for the scalar components, the scalar potential is then 

v( s , 0, ft = K 2 M 2 (M 2 + M 2 ) + ^ 2 |# - m 2 | 2 + ^ (e - \4>\ 2 + 1«^| 2 ) 2 , (4-ii) 

where g is the U(l) gauge coupling and a D-term with Fayet-Iliopoulos (FI) pa- 
rameter £ has been added. One can see that the charged fields become unstable at 
|s c | 2 = (M 2 + £<7 2 /k 2 ). If the FI term dominates the potential energy density, we 
have D-term inflation, 77 ) and if subdominant (or absent) F-term inflation (see e.g. 
78)). 
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The classical potential is completely flat at <f) = 4> = 0, with V = k 2 M 4 + \g 2 ^ 2 , 
and one needs to take into account the radiative corrections for the system to be 
able to evolve towards its minimum. For |s| 3> |s c | one has 

AV S = ^(^M 4 + ^ 2 )ln (^pj , (4-12) 

where fi is the renormalisation scale. In models where the inflation fields (5, 

are coupled to the supersymmetric standard model fields (see e.g. Refs. 83)) the 

potential is somewhat more complicated. 

This class of simple supersymmetric inflation models is highly predictive, as 
the cosmic string and inflationary scalar perturbation amplitude are to first order 
controlled by the same parameter M 2 +£<? 2 / k 2 . 79 ) One is pushed towards the F-term 
model with a very small coupling Kj 38 )' 80 ) Q r towards a non-minimal supergravity 
model. 38),82) 

4.3. Elementary cosmic strings 

Strings may be of any length, and so the question arises whether we could 
expect to find "elementary" strings of the size of the cosmological horizon, L ~ t. 
The question was considered by Witten. 84 ) Firstly, one can expect to find only 
the strings of a closed string theory: long open strings would fragment with a time 
of order the inverse string scale, M s _1 . Secondly, the strings will be very massive: 
Newton's constant is related to the string tension l/(27ra') by G ~ a', and so 
one expects Gfi = 0(1). Thirdly, cosmic superstrings would also be axionic, and 
therefore bounded by domain walls of surface density ^1q CD , where Aqqb ~ 1 GeV 
is the QCD scale. 

Following the discovery of D-branes, 85 ) large extra dimensions, 86 ) and warped 
compactification 87 ) the question was reconsidered. Importantly, the new compactifi- 
cations could reduce the effective 4-dimensional string tension well below the Planck 
scale. For example, in a warped compactification scenario the metric takes the form 

ds 2 = e 2A{y) n^dx^dx y + dyl, (4-13) 

where y a (a = 1 ... 6) are the compactified coordinates and (// = ... 3) are 
coordinates in order 4-dimensional spacetime. This means that a fundamental string 
living at y = yo, where the "warp factor" is Aq, has a tension 

"=w- < 4 ' 14 > 

and with sufficient warping the string tension can be well below the Planck scale. 

Strings from string theory come in two kinds: the fundamental or F-strings of 
the original theory, and D-strings, which are 1-dimensional D-branes, objects on 
which fundamental open strings end (see e.g. 88))*). They can form bound states 89 ) 



*^ To distinguish these strings from solitonic strings F-strings and D-strings are together called 
elementary in this review. 
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of p F-strings and q D-strings, with the composite (p, (/)-string having tension in a 
flat 10-dimensional spacetime 

fip,g = ^J"Jp 2+ 9s 2 Q 2 , (4-15) 

where g s is the string coupling, which is of order unity. More complicated formu- 
lae apply in the presence of background fluxes and in warped compactifications. 16 ) 
Provided these strings are well-separated in the extra dimensions from other, higher- 
dimensional, D-branes on which they can break, they are essentially stable. 

String theory therefore contains an interesting spectrum of objects which could 
constitute cosmic strings, and it also has a possible formation mechanism through 
the collision of D3 and anti-D3 branes at the end of brane inflation. 12 ) It turns out 
that this model looks rather like a hybrid inflation model, with the inflaton being 
the brane separation modulus, and the waterfall field the tachyon fields associated 
with the string states connecting the two branes. 

The low-energy effective action of a macroscopic quantum string is just the 
classical Nambu-Goto action with the string coupled to massless fields: perturba- 
tively, these are the spacetime metric, the dilaton, and a Ramond-Ramond scalar 
(equivalently a 2-form field). On phenomenological grounds, the scalars should gain 
a potential as part of the compactification process, and as a result become massive 
and fixed in value. Therefore we expect the strings with coordinate X^{a a ), a = 0, 1, 
moving in a 4-dimensional spacetime with metric g^ , to be described by the action 




(4-16) 



where ^y a p = g^J daX^daX 11 is the induced metric on the string worldsheet, and R 
is the 4-dimensional Ricci scalar. 

An important feature of elementary cosmic strings is the presence of junctions. 
A (p, q)-string can split into e.g. a {p — l,g)-string and an F-string at a three-way 
junction, a feature which has resulted in a lot of work on the classical dynamics of 
cosmic strings with junctions. 92 ) 

Certain field theories, including the Nielson-Olesen string at A/e 2 < 1, can 
also produce bound states and junctions, which means that the junction is not a 
distinguishing feature of an elementary cosmic string. On the other hand, these 
theories can be used to study the dynamics of networks of strings with junctions. 94 ) 

§5. Dynamics of cosmic strings 

5.1. Abelian Higgs model in an expanding background 

In Section 4.1 we exhibited the Nielsen-Olesen string solution of the Abelian 
Higgs model, a static and cylindrically symmetric field configuration representing 
a straight string. One can expect to find more general time-dependent solutions 
representing moving strings, either infinite or in the form of closed loops. If we wish 
to study these time-dependent solutions in a cosmological context, one must take 
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the action for the field theory in a Friedmann-Robertson- Walker (FRW) background 
metric, 

S = - j d^x^-g (g» v D^*D v <P + V{<j>) + ^g^g^F^F^ , 

where we have rescaled the gauge field — > A^/e. The background metric is 

ds 2 = a 2 {T)(-dr 2 + dx 2 ), (5-1) 

where dx 2 is the 3-dimensional line element, and a(r) oc r or r 2 depending on 
whether the universe is radiation or matter-dominated. 
In the temporal gauge (Aq = 0) the field equations are 

+ 2-^- D 2 cf> + \a 2 (\<j)\ 2 -<p 2 ) ( j) = 0, (5-2) 
(J^F^ ~ ia 2 (<p*D^ - D„if>*<f>) = 0, (5-3) 

where the index on the partial derivative is raised with the Minkowski metric. 

One can solve this system of partial differential equations numerically. 22 )' 23 )' 33 ) 
The scalar field is placed on a cubic grid with lattice spacing Ax, with the gauge field 
on the links. The boundary conditions are toroidal, and the lattice spacing must be 
chosen so that the core of the strings can be resolved: in practice Ax = 0.5m" 1 is 
sufficient. 

The initial conditions model an initial thermal or quantum vacuum state in 
which correlations are short-ranged: the simplest way to do this is to select the real 
and imaginary components of the scalar field from a gaussian distribution of zero 
mean and unit variance, independently at each site, at initial time t\ = Ax. The 
subsequent evolution generates a random gauge field, with the Hubble damping term 
removing the high frequency modes and cooling the system so that the scalar field 
settles into the minima of the potential at \4>\ = 4>q. 

There is a technical problem to solve: the physical width of the string is fixed by 
the mass scale m v , and so the comoving width of string shrinks, as a" 1 . Eventually 
the string will no longer be resolved by the lattice and becomes stuck on the grid. 
One can adjust the lattice spacing so that Ax = 0.5m" 1 at the end of the simulation, 
however, this means that the width of the string is large at the initial time and the 
simulation box can contain fewer strings. This is a practical strategy in the radiation 
era, where a oc r, but not in the matter era, where a oc t 2 . 

To get around the problem, one can modify field equations 33 )' 95 ) according to 

4> + 2-j> - D 2 <j> + Xa 2s (\<P\ 2 - 4)<P = 0, 

/ a 2(l-s) \ 

— ^— F„ v - ia 2 (4?D v 4> - D v <t)*<t>) = 0, 

If s < 1 the comoving width of the string does not shrink so fast, but the physical 
width of string w = m~ l "fattens" as a l ~ s . The string tension remains constant 
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as the ratio of the effective couplings is constant, which means that the collective 
dynamics of the strings is unaffected. 

This particular choice of modification preserves Gauss's Law (and hence cur- 
rent conservation) but because of the time-dependence violates covariant energy- 
momentum conservation. One must therefore carefully check the dependence of 
results on the adjustable parameter s, which as mentioned above can be taken to 
the physical value s = 1 in the radiation era. 




Fig. 3. Left: snapshot from a numerical simulation of a string network in the Abelian Higgs model, 
as described in the text. The lines show the centres of the strings, and the shading on the top and 
bottom faces represents magnetic field energy density and scalar field potential energy density 
respectively. Right: string network length scale £ as a function of conformal time, averaged over 
two sets of 5 simulations with different initial conditions (see Ref. 34) for an explanation). After 
an initial relaxation period the length scale evolves linearly, meaning that the network is scaling. 
The parameters of the simulation were: A = 1, e = 1, Ax = 0.5 at the end of the simulation, 
At — 0.1, volume 768 3 (upper line) and 1024 3 (lower line), with s = 0.3,0 respectively. 

A snapshot from a simulation on a 512 3 lattice in the radiation era is shown in 
Figure 3. The random initial conditions on the field have produced a random network 
of strings, most of which stretches across the box and so represents infinite strings. 
The shading on the top and bottom faces picks out the fluctuations in the gauge and 
scalar field energy densities, generated by the decays of the strings themselves. As 
the system evolves, the total length of string reduces, and all the energy eventually 
ends up in the form of gauge and scalar radiation. 

To be more quantitative, one can examine how the (comoving) length L of string 
evolves with time. It is convenient to look at the quantity 

S=y/iV/L), (5-4) 

where V is the comoving simulation volume, which can be thought of as the approx- 
imate distance between strings. Examining Figure 3, one sees that after an initial 
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relaxation period £ settles down to a linear growth in conformal time 

e oc r, (5-5) 

with the same constant of proportionality for two different ways of setting up the 
initial conditions. 34 ) One can measure the length in a number of different ways, but 
they are all proportional to each other. 

A linear growth of the comoving length scale with conformal time is easily seen 
to imply that the physical length scale £ p = a£ grows linearly with cosmological time 
t = f a{r)dr. Hence the physical length per unit volume is decreasing as t~ 2 . Given 
that the energy per unit length is fj,, the physical energy density in string goes as 

Ps ~ ^- (5-6) 

Hence the string density parameter Q s = 8TrGp s /3H 2 is constant, fi s ~ Gfi, a 
property known as scaling. 

5.2. Nambu-Goto strings 

In the limit that the ratio of the string width w to the curvature radius R goes 
to zero, one can treat solitonic strings as being ideal relativistic line objects, 96 )~") in 
which case they should be described by the Nambu-Goto action coupled to gravity 
(4-16), and also a massless pseudoscalar in the case of a global string. 62 ^ Large 
elementary strings are also described by this action. 90 ) 




Fig. 4. Worldsheet of an ideal string in spacetime, labelled by coordinates <r° (timelike) and a 1 
(spacelike). 

In studying the dynamics of Nambu-Goto strings it is convenient to chose world- 
sheet coordinates <r a (see Fig. 4) such that 

X-X = 0, X 2 + X 2 = 0, (5-7) 
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where the dot (prime) indicates a derivative with respect to a (er 1 ), and to choose 
the temporal or static gauge X° = a , identifying worldsheet time with 4-d time. In 
this gauge the equations of motion in a FRW spacetime become 

e + 2^(l-X 2 )e = 0, (5-9) 



/ ' 2 -2 

where e = y X /(l — X ). Here, e is the comoving energy per unit comoving length 
as measured by a = a , and can be chosen arbitrarily at some initial time. 

For short wavelength oscillations and on timescales much smaller than the Hub- 
ble time, the strings are effectively moving in Minkowski spacetime, for which we 
can choose e = 1, and the equations become 

X-X = 0, X-X = 0, X 2 + X 2 = l. (5-10) 

It is clear that the equations take the form of a wave equation, with unit propagation 
speed, supplemented by constraints which keep the velocity orthogonal to the tangent 
vector and the energy per unit length constant. 

When strings cross, the classical Nambu-Goto action must be supplemented with 
a rule determined from the underlying theory. Numerical simulations in the Abelian 
Higgs model 101 )' 102 ) show that crossing strings reconnect, unless they are moving 
very close to the speed of light, or are nearly parallel. Network simulations show 
that this is rather unlikely, and that the reconnection probability averaged over all 
angles and collision speeds is close to 1. Elementary strings on the other hand, being 
quantum mechanical objects, reconnect with a probability 

Pr~9 2 s ^-f(9,v), (5-11) 

where g s ~ 1 is the string coupling, f(9, v) depends on the orientation 9 and relative 
speed v, and V± is the volume of the extra dimensions in which they move. 103 ) This 
probability may be rather small, which is a key difference between elementary and 
solitonic strings: the lower the reconnection probability the higher the string length 
density. Numerical simulations in flat space are consistent with (L/V) oc P,." 1 , 104 ) 
while in cosmological backgrounds (L/V) ~ P~ ' 6 . 105 ) 




Fig. 5. Reconnection of a string leading to loop formation. Two pairs of oppositely-travelling kinks 
(discontinuities in the string tangent vector X') are formed. 
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When strings curl back on themselves and reconnect, they produce a loop, and 
creating a pair of oppositely-moving kinks (discontinuities in the string tangent vec- 
tor X') both on the loop and the parent string. Kinks are an important feature 
of string networks: classical Nambu-Goto network simulations 17 )' 18 )' 20 ) show that 
strings are full of kinks (see Fig. 6) left behind by loops. This loop production allows 
the infinite string length density to scale, although the total length of string is fixed 
by conservation of energy. 




t/i 



Fig. 6. Left: snapshot from a Nambu-Goto string simulation, showing loops formed during the 
relaxation from the initial state (blue) and during scaling (red). Right: infinite string comoving 
length scale £ (as computed from the invariant length L — f edo) for a Nambu-Goto simu- 
lation. 18 - ) Distance and conformal time are given in terms of the initial comoving correlation 
length l c , which is the only fixed length scale in the system (snapshot and data courtesy C. 
Ringeval) . 

Note that the slope of £ is lower than that for Abelian Higgs strings, meaning that 
the infinite string length density is higher, by about a factor of 2.5. 24 ) However, the 
main difference is clearly the presence of string loops in the Nambu-Goto simulations. 
The visual impression in Fig. 6 is dominated by the loops produced earliest, during 
the relaxation to scaling, which disguises fact that the loops are produced in such a 
way as to maintain a scaling power law distribution for loop size. 

Scaling arguments show that the number of (stable, non-self-intersecting) loops 
with physical length in the interval [1,1 + dl] is 25 ) 

n(l,t) = dj?n(l/t), n(x)cKx 3u ~ A , l>h = a(t)l c (5-12) 

where the scale factor grows as a oc t u , l c is the initial comoving correlation length, 
and dh is the horizon distance. *) 

One can define a loop production function by f(l,t) = a _3 Jj (a 3 n(/,t)) , and 
express it in scaling form 

f(l,t) = d^f(x), /(x)cxx-p. (5-13) 

*' Departures from the scaling power law indices were observed in, 18 ) which was suggested in 21 ' 
to arise from loop production at small scales, or to be a sign that the loop distribution was not yet 
scaling. 
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Recent simulations 21 ) are uncovering this scaling form, which takes a surprisingly 
long time to emerge, and show a very broad distribution consistent with p ~ 2 over 
a couple of orders of magnitude in x, with broad peaks at x ~ 0.05 and at the lower 
cutoff /* = a(t)l c set by the initial comoving correlation length l c . One can show that 
J xf(x)dx is bounded by conservation of energy, and hence p < 2. The true loop 
distribution must take into account the decay of loops into radiation, which modifies 
the simple power laws. 106 ) 

The difference in the loop number densities between Nambu-Goto and Abelian 
Higgs simulations demands some explanation. If one is treating the Nambu-Goto 
action as an approximation to the motion of solitonic strings, one is operating beyond 
its strict domain of validity (large curvature radius) as the string is full of kinks 
where the curvature radius R is of order the string width w. There is therefore 
no particular reason to expect the behaviour to be the same. In the Nambu-Goto 
approximation, the massive modes of the field are eliminated from the outset, and 
so infinite Nambu-Goto strings can lose energy only by producing loops. However, 
it is puzzling that Abelian Higgs strings with w/R <C 1 can produce radiation with 
frequency uj ~ w~ l ~ ecfro, while naively it appears to be made of fields changing with 
a frequency i? _1 <C w . Indeed, for a smooth string prepared in a standing wave 
of wavelength ~ R the power per unit length in massive radiation goes as dP/dL ~ 
(fj 1 /R)ex.p(—R/w), 12 ^ and massive radiation becomes negligible as R — > oo. Thus 
it is often argued that a network of Abelian Higgs strings must eventually start 
behaving like Nambu-Goto strings. However, real string networks do not prepare 
their initial conditions so carefully, and massive radiation continues at a rate per 
unit length dP/dL ~ (n/R) in the largest simulations to date, where w/R ~ 10~ 2 . 
In 2D domain wall networks massive radiation operates at this rate at least until 
w/R ~ 10" 3 , 107 - 1 so there is clearly a mechanism at work which cannot be ignored. 
While the mechanism is not well understood, small scale structure on infinite strings 
(discussed in the next section) has been shown to be important, 24 ) and such energy 
cascades from large to small scales are familiar from the study of turbulence. 

Thus it is an open question whether Nambu-Goto strings approximate the mo- 
tion of solitonic strings at late times, and in particular whether soliton strings pro- 
duce any significant population of loops. This has an important bearing on all 
observational signals from loops. 

For macroscopic elementary strings the situation is also unclear as the technology 
to do a proper quantum calculation is not well evolved: it is difficult enough even 
to construct the states. 108 )> 109 ) However, given that the low-energy approximation 
for a macroscopic elementary string is just the classical Nambu-Goto action, an 
intriguing possibility is that Nambu-Goto simulations are actually modelling cosmic 
superstrings rather than their solitonic counterpart. If this is true, it should be quite 
easy to tell solitonic and elementary strings apart based on observations of the loop 
population alone. 

5.3. Small scale structure 

Early Nambu-Goto string simulations 17 ^ showed that there were "wiggles" on 
infinite strings, excitations with wavelengths much less than the network length scale 



16 



Mark Hindmarsh 



£. The wiggles took the form of closely spaced kinks, with as much as 50% of the 
energy of the infinite string contained in these high frequency waves. Characterising 
the wiggles by their fractal dimension df on a scale I, numerical simulations show that 
df rises gradually from 1 on scales small compared with £ (approximately straight 
strings) to 2 on large scales (Brownian random walks). 19 ) 

This small scale structure can also be characterised by the tangent vector cor- 
relation function, which goes as 111 ) 

C T (a, t) = (x(a) • X(0)) = (±(0) 2 ) (l - A(a/t) 2 *) , (5-14) 

where A is an 0(1) constant and the exponent x determined from the mean square 
velocity v 2 and the Hubble expansion exponent v as 

This scaling power law form for small scale structure has been found both on Nambu- 
Goto 111 ) and Abelian Higgs strings. 24 ) At the very smallest scales, the correlation 
function approaches a constant as a, explicable in terms of the kinks. 112 ) The kink 
separation decreases very slowly, and scaling is absent. 

The significance of the small scale structure is that it is involved in transporting 
energy from large to small scales on infinite strings. The small scale structure triggers 
loop production near cusps (points on the string where the tangent vector vanishes 
and the string moves at the speed of light). 111 ) On Abelian Higgs strings it contains 
the necessary high frequencies for producing massive radiation. 24 ) 

5.4. Gravitational radiation reaction 

Early string simulations 17 ) showed no sign of the expected scaling in the loop 
distribution. The belief in scaling was strong, and it was conjectured that the typical 
loop size was set by gravitational radiation damping of small amplitude waves on 
infinite strings. 17 ) In the radiation reaction scenario the relevant scale for loop 
production is (I) ~ (G/i) p i, with p variously estimated to be l, 56 ' (1 + 2/3)/2, 114 ) 
and 1 + 2x, m ) where (3 = 1(2) in the matter (radiation) dominated era and x is the 
exponent of the power spectrum of small-scale structure on infinite strings. 

A numerical check of gravitational backreaction 115 ) demonstrated the expected 
damping of high frequencies on a simple kinky loop. However, no network simulation 
includes gravitational radiation reaction. In particular, it is not clear how damping 
operates on a scaling network which is self-intersecting and producing kinks. 

5.5. Summary 

There are two scenarios for ordinary cosmic strings without junctions, which we 
refer to as the Nambu-Goto and Abelian Higgs scenarios, to make clear their origin. 

Both recognise that the length density (both comoving and physical) in infinite 
strings scales, that is Lt 2 /V is a constant. In the Abelian Higgs scenario this is 
about 2.5 - 3 times as high as in the Nambu-Goto scenario (see Table, after 34)). 

In the AH model, strings lose their energy primarily into scalar and gauge ra- 
diation. There is 0(1) horizon-sized loop produced per Hubble time, and the loops 
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Simulation 


Measure 


Radiation 


Matter 


AH 34) 




0.255 ±0.018 


0.285 ±0.011 


AH 34 ) 


Lt 2 /V 


15 ±2 


12.2 ± 0.96 


NG 18J 


Lt 2 /V 


48 


27 


NG 19) 


Lt 2 /V 


37.8 ± 1.7 


28.4 ±0.9 


NG 21) 


Lt 2 /V 


44 


35 



Table I. Numerical results for the network length scale in horizon units and the string length density 
L/V normalised to the horizon size r. AH labels results from Abelian Higgs simulations, while 
NG labels results from Nambu-Goto simulations. Note the NG results quoted includes only 
infinite strings, and that not all authors include an error estimate. See 24 ' for a discussion of the 
AH string length measure. 

have a lifetime of order the Hubble time, decaying into radiation just as the infinite 
strings do. 

In the NG model, strings lose their energy into loops, whose sizes range between 
the comoving initial correlation length and about a tenth of the horizon size dh . The 
average loop size at formation is the comoving initial correlation length. On scales 
between the comoving initial correlation length and O.ldh the rate at which energy 
(length) is emitted into loops is approximately constant per unit logarithmic length 
interval. 

Gravitational radiation reaction is widely believed to play a critical role in setting 
the lower cut-off on a real string network. In view of the importance of establishing 
the typical loop size at formation for many observables, particularly gravitational 
radiation background from strings, a simulation including radiation reaction is re- 
quired. 

Finally, it should be noted that we are extrapolating from the end of inflation 
(U ~ 1CT 36 s) to today (i ~ 3 x 10 17 s). Over that range of scales even logarithmic 
corrections become O(100), so acquiring a better theoretical understanding of scaling 
in string networks, perhaps by techniques taken from other fields, 117 ) is extremely 
important. 

§6. Observational signals of cosmic strings 

6.1. Stochastic background of gravitational radiation 

An oscillating loop of string is the source of gravitational radiation with total 
power 54 )' 116 ) 

p = rGfi 2 , (6-i) 

where r ~ 10 2 is a dimensionless constant depending on the loop trajectory. The 
radiation power of a loop of length I is peaked at the fundamental frequency uj\ = 
4ir/l, and for high harmonics uj n decreases as n -4 / 3 , due to the strongly beamed 
radiation from cusps. Both the stochastic background from the integrated sum of 
all loops on our backward lightcone, 54 )' 58 )' 118 ) and the rare but strong signals from 
individual cusps 59 -* offer prospects for detection. 

Let us denote the average size at formation by l\ = at, and also to scale out a 
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factor G\i and define e = a/Gfi. In the original radiation reaction scenario e was a 
constant independent of G/i, but in more recent versions it is a higher power of G/i. 

There are limits on the stochastic background from millisecond pulsar timing, 118 ) 
from Big Bang Nucleosynthesis, 12 °) and from LIGO. 121 ) The strongest limits are 
from the pulsar timing, giving 1 19 ) 

Q g h 2 < 2 x 1(T 8 (6-2) 
In the Nambu-Goto model, the stochastic background from string loops is, for a <C 

where x* = £ p /dh and J? m is the matter density parameter today. Hence one can 
estimate G[i < 1CP 6 . A more accurate analysis 41 ) gives 

Gp<7 x 10~ 7 (6-4) 

As there is only 0(1) loop per Hubble volume in the Abelian Higgs model, there are 
no constraints from gravitational radiation. 

6.2. Cosmic rays from cosmic strings 

Solitonic strings can decay into the particles of the fields from which they are 
made. There are three main avenues: perturbative production; 122 )' 123 ) cusp annihi- 
lation; 124 )' 125 ) and non-perturbative massive radiation. 22 )' 126 ) 

For loops found in network simulations of size I and width w ~ m^ 1 , the power 
in the three mechanisms scales as 

1. Perturbative particle production: P ~ fj,(w/l) 

2. Cusp annihilation: P ~ ^(w/l) 1 ^ 2 

3. Massive radiation: P ~ jjl. 

The largest flux, and therefore the most important bounds, comes from non- 
perturbative massive radiation, assuming that the scaling configuration observed in 
AH simulations continues indefinitely. Cosmic strings produce massive particles, 
denoted X, at a rate 

n x ^ — , (6-5) 

m x \t J 

where to is the time today, and the energy injection rate Qo is given by 



6// 

x 2 d h (t ) 



Qo = -Ps(t ) = —j u , 3 . (6-6) 



Cosmic strings constitute a p = 1 topological defect (TD) 66 ) model, for which there 
is an upper bound on the energy injection rate from the low energy diffuse 7-ray 
background 66 ) 

Qo < 4.4 x 10~ 23 /i eV cm" 3 s -1 . (6-7) 
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Assuming the massive particles decay into Standard Model particles with a non-zero 
branching fraction /, one can derive a bound 70 -*' 72 )*- ) 

G/m < lO' 10 /' 1 . (6-8) 

Although the strings produce particles with energies up to mx, which may be up to 
10 16 GeV, there is no bound from Ultra-High Energy Cosmic Rays (UHECRs, defined 
as cosmic rays with energies greater than 10 20 eV). The range of these particles is 
less than 100 Mpc, and the nearest string in the solitonic string scenario of order a 
Hubble distance**^ 

If strings are described by the Nambu-Goto model, there is no significant bound 
from cosmic rays. Hence cosmic ray and gravitational radiation bounds are comple- 
mentary. 

6.3. Cosmic Microwave Background constraints 

The gravitational field of a moving string between us and the last scattering sur- 
face produces a discontinuity in the apparent temperature of the Cosmic Microwave 
Background (CMB), 26 ) 

AT ~ 87r(Gfi)vT C MB, (6-9) 

where v is the transverse speed (see Figure. 11). Hence (AT /T) ~ 10~ 5 (G///10 -6 ), 
which means that GUT or inflation-scale strings produce observationally interesting 
CMB perturbations. To resolve the distinctive edges of the Gott-Kaiser-Stebbins 
(GKS) effect one needs both high resolution and high sensitivity. A number of limits 
can be derived by directly looking for these edges: from OVRO 127 ) Gfj, < 11 x 10 _6 x*, 
where x* ~ 0.3 is the string length scale in units of the horizon, and from WMAP 128 ) 
G\i < 3.7 x 10~ 6 . See also Refs. 129), 130) for other edge detection methods. 

More stringent constraints come from calculations of the CMB temperature and 
polarisation power spectra and their comparison to the data, principally WMAP. 33 )' 37 )' 38 )' 41 ) 
To calculate power spectra, one solves the Boltzmann equation for the photon-baryon 
fluid coupled to the linearised Einstein equations in an expanding background, with 
an additional source in the form of the energy-momentum of the cosmic strings 
(see e.g. Ref. 132)). Schematically, the Fourier-space perturbation equations can be 
written 

V a fs(T,k)h^T,k) = S a (T,k), (6-10) 

where h a (r,k) are the perturbations in the metric, matter, radiation intensity, or 
polarisation; D a p(r,k) is a (time dependent) time evolution operator; and S a (r,k) 
is the source. These equations can be re-expressed as coupled evolution equations 
for the multipole moments of the CMB angular power spectrum. In the standard 
cosmological model, inflation determines the initial conditions for hp(r, k) and there 
are no sources. Standard Einstein-Boltzmann solvers 28 ) evolve the initial conditions 
forward and deliver the CMB angular power spectra, and also the matter power 
spectrum. 

*' Note that these authors give a bound a factor of 10 lower, as they took a value of the string 
density approximately 10 times higher than indicated in current Abelian Higgs simulations. 
**' The bound derived in Ref. 22) from the UHECR flux is therefore incorrect. 
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The CMB angular power spectra are derived from the multipole moments of 
the temperature and polarisation distributions. For example, for the temperature in 
direction n one defines a\ m = f df2AT(n)Y,*(n), and the angular power spectrum 
is 

l 

Cl=Yl \ a lm?- (6-11) 

m= — l 

Most often plotted is the anisotropy power, 1(1 + l)C//(2-7r). On a linear-log plot the 
power over a given multipole range is proportional the area under the curve. 

The standard techniques can be adapted to include sources. The contribution 
to the power spectra from the sources can be written schematically as 134 )> 135 ) 

(\K(r ,k)\ 2 ) = j j V 1 ©" 1 (S a (T,k)S* a (T> \k)) (6-12) 

where the angle brackets indicate an average over an ensemble of source histories, and 
we denote the numerical integration symbolically by T> . There are two approaches 
to computing the averages. 

In the unequal-time correlator (UETC) method one performs a number of nu- 
merical simulations, and averages to find the UETCs of the source energy momentum 
tensor 

C^ p x(k,T,T') = (T^(k,T)T* pX (k,T')) . 

Due to the rotational symmetry of the system, there are in fact only 5 independent 
UETCS: 3 scalar, 1 vector, and 1 tensor. Furthermore, scaling makes the correlation 
functions behave as 

z n 

C(k,T,r') ~ 4—=C(z,r), (6-13) 

V TT 

where r = r/r', z = k\J tt', and n = 2 is conventionally chosen for the vector 
correlator. The dimensionless correlators C are peaked at r = 1, z ~ 1 with a ridge 
at r = 1 of width 0(1) for z < 1, and 1/z for z > 1. In practice scaling is vital 
to making the calculation feasible: simulations over only a few expansion times are 
required. Figure 7 shows the vector UETC from an Abelian Higgs cosmic string 
simulation. 34 ) 

The UETCs are symmetric under interchange of z and z' , and therefore can be 
decomposed into orthogonal eigenvectors ^\(k,r) of the UETCs, defined from 

° dx'C{x,x')ip x (x') = \^J X {x), (6-14) 

where x = kr, xq is the value of x today, and i; < k One can then apply the 
Einstein-Boltzmann solver to the eigenvectors, and recover the power spectra by 
summing over the separate contributions. 

Another approach is to model the sources more economically than a full 3D 
simulation of the field equations, in such a way that the UETCs are reproduced. For 
strings, the Unconnected Segment Model 36 )' 67 ) models the strings as moving sticks 
of energy fx, tension T and length L(t) = x*t. The sticks have random positions and 
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Fig. 7. The vector unequal time correlator (UETC) from a cosmic string simulation. 34 ' Light grey 
comes from an average over data from three simulations, dark grey indicates extrapolation. 



velocities v, and decay at random times so that the energy density decays as t~ 2 . 
It has four parameters: fi, the ratio f3 = T/fi, and v, which can be related by 
a phenomenogical model such as the velocity-dependent one-scale (VOS) model. 137 ) 
There is enough freedom to model either the Abelian Higgs or Nambu-Goto simu- 
lations 41 ) (see Fig. 8). Nambu-Goto simulations are modelled with larger /3 and v, 
and smaller a;*, and tend to have a relatively smaller vector contribution. 
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Inflation+strings 
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Fig. 8. Left top: Strings normalised to WMAP 3 {£ = 10) 33) Left bottom: Differences from best-fit 
standard ylCDM cosmology. Right: Comparison of scalar (red), vector (blue) and tensor (green) 
power spectra in the Unconnected Segment Model (solid) with those produced from the UETCs 
of the Abelian Higgs model simulations (dotted). 41 ' 



Armed with the power spectra, one can apply standard Markov Chain Monte 
Carlo (MCMC) techniques (see e.g. Ref. 138)) to fit models with both inflation and 
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cosmic strings to the data. 33 ) The two sources of perturbations are not correlated: 
they are produced by fluctuations in different fields at different times, and so the 
power spectra can be simply added. The cosmic string power spectra are proportional 
to (G/i) 2 , so this results in its simplest form in a cosmological model with one 
more parameter, adding to the standard set f2bh 2 (fractional baryon density), Q c h 2 
(fractional Cold Dark Matter density), 9a (acoustic scale), V(ko) (scalar amplitude 
at the pivot scale ko), and n s {ko) (scalar tilt). One can also usefully parametrise the 
string contribution by its fractional contribution to the temperature angular power 
spectrum at a multipole I = 10, /io = C^ T ' string /C^ T ' total . 

One finds a degeneracy between G/j, and n s (or Qbh 2 ), and it was found that the 
best fit to the CMB data (3-year WMAP, with further CMB data on small scales) was 
a mixed inflation and AH string model with /io = 0.099 with n s = 1.00. 33 ) The de- 
generacy can be lifted by fitting to data such as the SDSS luminous red galaxy (LRG) 
power spectrum, 45 ^ or applying priors from Hubble parameter measurements, 139 ) 
Baryon Acoustic Oscillations, 46 ) and deuterium abundance measurements. 14 °) The 
CMB remains the most robust from both the observational and theoretical point of 
view. 

The most recent results from fits to the CMB alone give G/Jt < 6.8 x 10~ 7 for 
the Abelian Higgs model and G\x < 2.8 x 10~ 7 for the Nambu-Goto model, 41 ) using 
the 5-year WMAP data release in combination with CMB data on smaller scales. 
Combining the 7-year WMAP data with recent ACT data gives Gp, < 1.6 x IO" 7 for 
the NG model. 43 ) 

In hybrid inflation models the string tension and the scalar amplitude are related 
by the parameters in the potential, which can therefore be directly constrained. In 
the minimal D-term and F-term inflation models the important parameters are the 
scalar coupling k and the symmetry-breaking scale M (see Eq. (4-10)). Figure 9 
shows the la and 2a likelihood contours derived from a MCMC fit to CMB data 
(including 5-year WMAP) and SDSS LRGs. 38 ) The best fit models give Q h h 2 in the 
range 0.024 - 0.025, which is in tension with the value inferred from the deuterium 
abundance D/H in damped Lyman-a systems, Q\Jp? = 0.0213 ± 0.0010 via light 
element abundances predicted by standard Big Bang Nucleosynthesis (BBN). 140 ) 
However, given that there are only 7 measurements, with abundances corresponding 
to fi^h 2 in the range 0.017 to 0.030, the BBN constraint should be treated with 
caution. 

Future CMB data will be sensitive to lower values of G[i. Planck is scheduled 
report its power spectra in early 2013, and there are plans for a dedicated satellite 
missions to measure accurately the B-mode of the polarisation power spectra. 141 )' 142 ) 

An important question is whether a string contribution can be distinguished from 
other departures from the standard adiabatic scenario, such as gravitational waves 
(predicted in large field inflation models), and other kinds of topological defects such 
as textures and semilocal strings. 143 )' 144 ) 

It is conventional to parametrise the gravitational wave contribution in term 
of the tensor-to-scalar ratio r = Vt/Vt- From Fig. (10) one can infer that Planck 
can distinguish a string model fio > 0.02 from inflationary tensor fluctuations. 143 ) 
The reach in /io for a future B-mode satellite depends on assumptions made about 
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Fig. 9. Marginalised likelihoods for D-Term (red, smaller values of k) and F-term (blue, larger 
values of k) models of Section 4.2. The string models are for Nambu-Goto (left pair) and 
Abelian Higgs (right pair) using CMB and SDSS data (from the second of Refs. 38)). 

foregrounds and the efficiency of separating out the contribution from the lensed 
E-mode power spectrum. With conservative assumptions, the B-mode missions will 
be sensitive to /io ^ 2 x 10~ 3 . 144 ) Under the assumption of perfect delensing, one 
can reach f 10 > icr 5 . 145 )> 146 ) 




Fig. 10. Likelihood contours (68% and 95%) of the marginalized 2D posterior distribution for a 
cosmological model with (a) tensor-to-scalar ratio r — 0.04 and no string component, (b) r = 
0.04 and defect fraction /io = 0.008, and (c)/io = 0.01 with no tensor component. 143 ^ The 
actual models are depicted with a cross. 

6.4. Cosmic string CMB non-Gaussianity 

The GKS effect introduces a highly non-Gaussian signal at small angular scales, 
where the string signal is expected to dominate the power spectrum. 

Let = AT/T. In the flat sky approximation, and in the light cone gauge 
X° + X 3 = t, one can show 127 )' 147 ) 

-k 2 O k = iSnGfik* [ daX A (a)e ik ' X ( a \ 



where A = 1, 2 are the 2D transverse coordinates. 

By taking an average over an ensemble of strings, either analytically or nu- 
merically, one can show that the resulting power spectrum has an approximately 
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i^ 1 tail, 27 )' 127 )' 148 ' 1 ' 149 ) which dominates the string-induced acoustic oscillations for 
£ > 3000 or angular scales of less than about 3'. 34 )> 40 ) However, calculations which 
include the acoustic oscillations with an Einstein-Boltzmann solver 35 ) show that they 
may tend to obscure the GKS edges (see Fig. 11). 




Fig. 11. Small angular scale maps of the CMB fluctuations induced by cosmic strings. Left: GKS 
effect only, size 7.2°, resolution 0.42". 148 ' Right: acoustic oscillations included, size 3°, resolu- 
tion resolution 0.70'. 35) 

One can also work out the CMB bi- and trispectrum at small angular scales from 
the GKS effect by averaging over string ensembles. 150 )' 152 ** The bispectrum bk lt k 2 ,k 3 
and the trispectrum T klk2k3ki are defined by 

(0 kl 9 k2 e k3 ) = 6 fcl , fe2 , fc;! (27r) 2 <5(fci + k 2 + k 3 ), (6-15) 
(6> fcl 6>fc 2 6>fc 3 6'fc 4 ) = T fclife2jfc3>fc4 (27r) 2 <5(fci + k 2 + k 3 + k A ), (6-16) 

and one finds that the GKS effect gives an equilaterial bispectrum of order 

b kkk ~ -(4 x 10- 14 )c ^^^ ( G/i - V. (6-17) 
v ; V A £2£6 ^0.7 x 10~ 6 J y ' 

where v, t, L£/A, and cq are dimensionless constants characterising the string net- 
work. All but Co are 0(1), with cq ~ 10 _1 . *' The GKS bispectrum (and therefore 
the skewness) is negative, and of order C^ 2 . It is not straightforward or unambiguous 
translating the predictions for the bispectrum into an equivalent local non-linearity 
parameter f n i. 153 ^ Using the relation bin — (2 x 10 _17 )/~ 4 / n i, 153 ^ one finds that 
\f n i\ ~ 10 3 co(G/x/0.7 x 10~ 6 ) 3 at ~ 1 or t ~ 1000. A comparison using a sum 
over multipoles gives ~ 75c (G///0.7 x lO -6 ) 3 . 152 - 1 

*-' The smallness of Co can be understood as the result of approximate time reversal invariance 
for string motion on small scales (B. Wandelt, private communication, see also Ref. 151)). 
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The trispectrum goes as 



(6-18) 



with p = 6(l + 2x/3), which is of order Cf. An interesting feature of the trispectrum 
is that the index p depends on the small-scale structure parameter x, and is therefore 
dependent on the mean square velocity of the string. Observational limits are not 
yet strong enough to provide competitive constraints. 

6.5. Strings and 21cm radiation background 

Neutral hydrogen between us and the last scattering surface, between redshifts 
30 < z < 300, absorbs CMB radiation at 21cm in the local rest frame, and the mea- 
surements of the radio background therefore offer a sensitive probe of the distribution 
of the matter in this range of £. 154 )~ 157 ) Cosmic strings contribute to the power spec- 
trum of 21cm fluctuations, 158 ^ although at a rather low level compared to the signal 
from standard adiabatic perturbations (see Fig. 12). The larger the collecting area 
the lower the Gp which can be distinguished: with over 100 km 2 (and a precise 
knowledge of the other cosmological parameters) one can reach Gp ~ 10~ 10 . 158 ) 







WMAP/Pulsar limit 






Limit predicted by ^ — — — ■ 
D-brane inflation — — 




21 cm full sky probe limit 







Fig. 12. Left: angular power spectra from inflationary adiabatic initial conditions (COBE normal- 
ized) and cosmic strings. Right: constraints from current and future experiments on Gfj, for a 
sky fraction of 10%, bandwidth of 0.4 MHz and integration time of 3 years. 158 ' 

String- induced 21cm background perturbations are cross-correlated with the 
CMB, 159 ) but the prospects for detection even with future instruments like the 
Square Kilometre Array (SKA) look weak. 

As in the case of the CMB, strings produce distinctive real-space features in 
the 21cm sky. 160 -* As they move through the neutral matter, a shocked wake of hot 
hydrogen accumulates behind them, leaving wedge-shaped regions of higher emission. 
Wakes at z ~ 30 from strings with Gp ~ 6 x 10~ 7 should be observable with SKA. 160 ) 



§7. Strings - the next generation 



We have seen that there are well-motivated inflation models which produce cos- 
mic strings. There are many distinctive signals from such models, over and above the 
standard inflationary ones. Reliable calculations require a quantitative understand- 
ing of the evolution of a network of cosmic strings from the end of inflation until 
today. There are two principal routes: firstly with simulations of a classical field 
theory such as the Abelian Higgs model, which contains all the important non-linear 



26 



Mark Hindmarsh 



physics of the string width, and secondly with the classical Nambu-Goto action, 
with which significantly larger volumes and longer times can be simulated. The two 
models give quantitatively different results for the infinite string density and diverge 
considerably on the loop number density and size distribution. 

The best understood signal is the CMB, which is also the best measured. De- 
pending on which model one uses to simulate the string network, and which datasets 
are used, current observations give 2a upper bounds on the string tension parameter 
G\i in the range (1.6 - 6.8) x 1(T 7 (see Section 6.3). Planck will improve sensitivity 
to G[i by a factor of 3 or so, and future B-mode missions such as CMBpol or COrE 
by a further factor of 3. 

Strings are also well constrained by either the limits on the stochastic grav- 
itational radiation background from millisecond pulsar timing (in the traditional 
Nambu-Goto scenario, Section 6.1) or by the soft 7-ray background (in the case of 
Abelian Higgs strings, Section 6.2). 

An urgent task for the near future is to include gravitational radiation into 
network simulations, to check the assumption in the standard scenario that the loop 
size at formation is set by radiation damping. The amplitude and spectrum of the 
gravitational radiation both strongly depend on this quantity and current fits to the 
data leave it as a free parameter, greatly reducing the constraining power. 

It is also to important to understand whether solitonic strings continue to radi- 
ate massive particles as the field theory simulations show, or whether there is some 
critical network length scale at which they start to behave like Nambu-Goto strings. 
Nambu-Goto simulations have made impressive gains in dynamic range and the ex- 
pected scaling is finally being seen loop production; however, the NG approximation 
assumes from the outset that massive radiation is unimportant for string networks, 
and so scaling behaviour in NG simulations does not in itself prove that massive 
radiation from solitonic strings will eventually turn off. One way to resolve the issue 
would be to observe NG-like behaviour in a sufficiently large AH string simulation. 
However in current simulations classical NG string networks do behave differently 
from classical AH ones. 

Large-scale structure from strings is also an area for development. String matter 
perturbations feature wakes, which are non-linear as soon as they are created, and 
therefore the growth of structure should be investigated using N-body simulations. 

Cosmic strings are features of several models which accommodate both current 
particle physics and cosmological data, and there are non-trivial constraints from 
the combined datasets. Cosmic string phenomenology is therefore part of the wider 
effort to map out a theory which accounts for all data. 

It is notable that the most significant difference between elementary and soli- 
tonic strings, the small reconnection probability, was discovered as the result of a 
quantum calculation in string theory, which strongly motivates the drive for a more 
string-theoretic approach to cosmic strings. String theory can give the decay rate 
of loops of string into both massless and massive states, including backreaction, 161 ) 
and therefore one can check the standard assumptions using states corresponding to 
realistic cosmic string configurations. 109 ) Calculations for networks of infinite F and 
D-strings remain a project for the future. 
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